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For entangled states of light both the amount of entanglement and the sensitivity to noise gen- 
erally increase with the number of photons in the state. The entanglement-sensitivity tradeoff is 
investigated for a particular set of states, multi-dimensional entangled coherent states. Those states 
possess an arbitrarily large amount of entanglement E provided the number of photons is at least 
of order 2 2E . We calculate how fast that entanglement decays due to photon absorption losses and 
how much entanglement is left. We find that for very small losses the amount of entanglement lost 
is equal to 2/log(2) » 2.89 ebits per absorbed photon, irrespective of the amount of pure-state 
entanglement E one started with. In contrast, for larger losses it tends to be the remaining amount 
of entanglement that is independent of E. This may provide a useful strategy for creating states 
with a fixed amount of entanglement. 



I. INTRODUCTION 

Entanglement lies at the heart of many quantum in- 
formation processing protocols Q. In order to generate 
entangled states one typically envisages a process that in 
the ideal case would lead to the desired state with the 
desired amount of entanglement. Imperfections such as 
noise and losses, however, degrade the quality of the state 
and diminish the amount of entanglement created. As a 
consequence, the task for which the entangled state is 
created (say, teleportation 0) will be accomplished only 
partially. An alternative strategy for generating a cer- 
tain amount of entanglement would be the following: use 
a process that in the ideal case would generate (much) 
more entanglement than needed, but that produces the 
desired amount in the realistic imperfect case. 

If it were increasingly difficult to generate larger 
amounts of entanglement, that strategy would not be 
very useful. Fortunately, it turns out that for light beams 
a Standard optical nonlinearity [2(| allows one to gen- 
erate more entanglement when the interaction time (or 
length) r is shorter. one generates log 2 M ebits in a 
time r = 1/M |5j. This effect exploits some of the pe- 
culiar properties of entanglement in infinitc-dimensional 
Hilbert spaces [(|. We may use that effect in the follow- 
ing way: take a very short time (or length) r = 1/M to 
generate a bipartite MxM dimensional pure state of two 
modes with log 2 M ebits of entanglement, with M large. 
Since the interaction time is assumed short we may ne- 
glect decoherence during that part of the procedure and 
consequently the state at that point is still pure. By the 
time we will actually use the state for some quantum- 
information processing protocol, there will have been de- 
coherence, turning the pure state into a mixed state p. 
The most relevant type of decoherence for light beams is 
that due to photon absorption losses, and that is the only 
type of decoherence we will consider here. We are inter- 
ested in the amount of bipartite entanglement between 
the two modes @ in the state p. 

Since p will typically still be a very high-dimensional 
state, we use a measure of entanglement that can actually 



be computed efficiently even in that case: the logarith- 
mic negativity En (p) y| . That quantity can be linked to 
two operational measures of entanglement, demonstrat- 
ing the possible use of the mixed-state entanglement in 
p. On the one hand Ej^(p) determines a bound (through 
the negativity in fact) on the fidelity that can be achieved 
in teleportation using a single copy of p, on the other it 
gives a bound on the amount of distillable (pure-state) 
9] entanglement present in many copies of p (see Ref. [8j 
for details). 

An alternative method to create certain desired entan- 
gled states of light was discussed in [ïo|: that method 
is very similar in spirit to that analyzed here but still 
relies on producing a pure state. The first part of the 
procedure is the same as above and consists of generat- 
ing an M x M-dimensional pure entangled state (with 
M large) of two modes. But subsequently one performs 
appropriate local measurements on both modes to (ap- 
proximately and probabilistically) project the state down 
to a 2 x 2-dimensional state of the form 

\a) ® \a) +exp(i<f))\ - a) ® | - a), (1) 

where \a) is a coherent state, and where the amplitude 
a is not too small, such that \a) and | — a) are to a very 
good approximation orthogonal. In that case, the states 
(JTJ possess one ebit of entanglement for any value of the 
)hase c/> and such states can then be used for teleportation 
n|, quantum computing I12J. quantum communication 
13| . and testing nonlocality [14|. 

Generally speaking, the sensitivity to noise increases 
with the number of photons in the state. This was noted, 
in particular, in experiments on generatin g Sc hròdinger 
cat states of the form \a) + exp(i<^)| — a) |l5l . lïii. It is 
also true for the related entangled coherent states |lll·IÏ7| . 
and this issue was subsequently investigated in more de- 
tail and generality for light beams in [ïg. For an even 
more general discussion of decoherence of macroscopic 
and microscopic objeets, see for example |l9j |. Thus, 
while increasing M leads to more entanglement it also 
leads to the state being more sensitive to decoherence, as 
the number of photons needed to create that much entan- 
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glement increases as M 2 . That is, one expects the state 
to lose entanglement at a rate that increases with M. 
Thus, one may expect an optimum number of photons 
for fixed interaction time and fixed amount of decoher- 
ence to exist that maximizes the amount of entanglement 
left after decoherence. Confirming that expectation is 
another aim of this paper. 



II. ANALYTICAL RESULTS 

A. Multi-dimensional entangled coherent states 

We hrst review some of the results of [5J that we need 
here. Suppose one starts out with a coherent state \y2a) 
and one lets the state propagate through a Kerr non- 
linear medium, described by a Hamiltonian of the form 
H = a) 2 a 2 for a time r equal to 1/M, with M an in- 
teger. Subsequently one splits the resulting state on a 
50/50 beamsplitter with the vacuum. The state of the 
two output modes so generated is of the form 



I*) = 



M 



(2) 



where the states are coherent states with amplitude 
a and phases —2nq/M, 



= \ae-xp(-2iriq/M)}. 



(3) 



The states © generalize the states Q and may be called 
multi-dimensional entangled coherent states. For ease 
of notation, here and in the following we suppress the 
dependence of the states and similar states defined 
below on the parameters a and M. Up to irrelevant 
ovcrall phase factors the coefficients f q are given by 



u = 
h = 



1 exp( ^+l) ); Modd (4) 



M 



exp (----), M( 

M M ' 



(5) 



The states © potentially possess log 2 M ebits of entan- 
glement, namely when a is sufhciently large to make the 
states |aexp(— 2niq/M)) and \a exp(— 2niq'/M)) nearly 
orthogonal for q ^ q' . This requires that |ct| 2 be at least 
of order M 2 . More precisely, if we choose a small number 
í< 1 the requiremcnt 



1(^1* 



9+1/ 



= 6, 



gives 



4tt 2 



M 



(6) 



(7) 



for large M. This follows from the expression for the 
overlap of two coherent states with complex amplitudes 
a and /3 



ce, 



exp(a*/3 - 



2 /2- 



72). 



(8) 



Thus, the shorter the interaction time r = l/M, the 
more entanglement the states © contain, provided the 
number of photons is sufhciently large. 



B. Decoherence due to photon absorption 

The decoherence of a mode due to photon absorption 
losses can be described by a parameter r/ that gives the 
fraction of photons surviving the absorption process. The 
loss process is determined by its effect on coherent states 
by 



\a) <g> |0> B -> \y/rja) <g> \yJ\-r}d) E , 



(9) 



where E refers to the environment, which will be traced 
out. We have a bipartite state on two modes, and when 
we assume both modes to decohere in the same way, the 
pure state |$) is turned into a mixture of the form 

M M 
9=1 P=l 



1*9) (*. 



V\ 



(10) 



Here we defined two new states, related to those defined 
in (J3J|. First, due to decoherence the amplitude a is re- 
duced by a factor of 77, and hence we used the states 



\è q ) = \^jacxp{^2TTÍq/M)), 



(11) 



to expand p in. Second, the state of the environment can 
be expanded, likewise, in terms of the states 



l*9> = \y/l -vaexp(-27riq/M)). 



(12) 



After tracing out the environment the overlaps of such 
states appear in the expression for p, hence the appear- 
ance of the factor ((^pl*,)) 2 in (JTDJ|. 



C. Decoherence of entanglement 

In order to determine the entanglement in the state p 
as a function of 77, ce, M, we use, as announced, the log- 
arithmic negativity as measure of entanglement 0- We 
first have to express the density matrix in some orthog- 
onal basis. This can be accomplished as follows: first 
definc the Al-by-M Grani matrix G by its entries 



G Pq = (* P I*9 



(13) 



Then suppose we have chosen M orthogonal states \xí 
for i = 1 ... M in terms of which we expand 



M 



(14) 



Similarly, we expand the density matrix p in terms of the 
basis formed by the set 



\xi)(xj\ ® |a; fe )(x;| 



(15) 
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The matrix A can be expressed in terms of the Gram 
matrix by using the relation 



AI 



Gpq — Aqj Ajp ■ 

Namely, in brief notation, we have 



.4 



í G f , 



(16) 



(17) 



with the superscript T denoting the transpose. Here we 
used the fact that G is a Hcrmitian matrix with positive 
eigenvalues, so that A, too, can be chosen as a Hermitian 
matrix with positive eigenvalues. Defining a matrix C by 
its elements 



C, 



mn,q 



fqAq rn Aq n 



we can write the density matrix p as 

p = CG 2 C T , 
where G2 is the matrix defined by 

(G 2 ) Pq = ((%\* q )ï 



(18) 



(19) 



(20) 



The partial transpose of p with respect to the first mode 
and with respect to the basis (|15|) is denoted by p Tl and 
is defined through 



Pkl .ran 



prnl.kn • 



(21) 



The logarithmic negativity of p, finally, is given by 

E N (p)=log 2 \\p T % (22) 

with ||.|| denoting the trace norm, which for Hermitian 
operators is equal to the sum of the absolute vàlues of 
the eigenvalues. The result Il22[l is exact and can be cal- 
culated numerically in a straightforward way by using 
Standard linear àlgebra. On the other hand, the size of 
the matrix p may become too large for large M to be use- 
ful in numerical computations. Here we use Matlab, and 
whenever M < 50 we use the exact result (|22|l . In other 
cases we can use certain approximations vàlid under the 
appropriate conditions. 

D. Some approximations 

In the case that the states |<E> g ) are almost orthogonal, 
we can approximate En as follows. First, the partial 
transpose of p is approximated by 

M M 

9=1 p=i 



|èp>($,|®|$,><$ 



pr 



(23) 



Then the M eigenvalues of this matrix can be written 
as 

AÍ_ = ±^exp(-(l-ry)|a| 2 |exp(2^(fc-TO)/M)-l| 2 ) 



for (1 < k < m < M) 
\ k = — for (1< k < M). 



(24) 



Hence 



where 



E N (p)^\og 2 [l + F] 



(25) 



M-l 



F = exp(-(l - ï])\a\ 2 \exp(2mk/M) - 1| 2 ). (26) 



fc=i 



We can further approximate this when 
e = 1 — 77 



(27) 



is a small number, e <C 1. In fact, when not only e is 
small but e\a\ 2 is small as well, we can simply expand 
the exponentials appearing in the expression for F to 
obtain 



E N {p) » log 2 [M-2e|a| 2 (AÍ-l)] 



log 2 M- e\af 



M-l 



log(2) M 



(28) 



When e becomes larger, such that e\a\ 2 is no longer small, 
we can make a different approximation in the case that 
relation holds with 8 <C 1 some fixed small number. 
In this case, F can be rewritten as 



M-l 



F=Y^ cxp(-elog(l/<5)M 2 |exp(27nífc/M)- 1| 2 /(2tt) 2 



fe=i 



(29) 

This sum cannot be performed analytically, but we can 
approximate it as 



M-l 



F 



fe=i 



(30) 



which is vàlid for e not too small. Since S -C 1 only 

terms with sufliciently small vàlues of k contribute. In 

particular, if we simply cut off the summation when 
,2 

5 e becomes smaller than some number D, i.e., at k = 
y/log(D) I log((5) y/ï/e, then the entanglement En{p) is 
of the form 



E N {p) « log 2 [l + yjD>/{e\a E {5))\. 



(31) 



We checked this approximation numerically and D' was 
determined to be around D' « 7T. 
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FIG. 1: Loss of entanglement per photon cLEn/cLN as a 
function of log 10 e. Two curves on the right correspond to 
M = 200, the two on the left to M = 2000. The sòlid curves 
correspond to 5 = 10~ 2 and the dashed curves to 5 = 10~ 4 . 
This plot uses the approximation 1251 . 
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FIG. 2: Change in logarithmic negativity, AEn as a function 
of the change in number of photons, AiV. There are in fact 
5 curves here, several of which overlap to make them indis- 
tinguishable. The dashed curve is for M — 20 and 8 = 0.2. 
The sòlid, almost overlapping curves, are for 8 — 0.01 and 
M = 20, and for 8 = 10~ 4 with M = 20, 200, and 2000. 



III. NUMERICAL RESULTS 

Here we use the exact result (12211 and the approximate 
results <|25|) . (|28|) . and (|31|) to investigate the behavior 
of entanglement as functions of various different param- 
eters. 

When a is sufficicntly large and r\ = 1 one will have 
log 2 M ebits of entanglement. But how quickly does one 
lose that possibly large amount of entanglement when 
there is a small amount of decoherence? The average 
number of photons absorbed from each mode is 



AiV = e\a\ 2 . 



(32) 



The change in entanglement from pure to mixed state, 



AE N = E N (p) - E N (\$)($\), 



(33) 



is, for very small losses, proportional to AiV and is given 
by 2/log(2) w 2.89 ebits per photon: 



AE 



M — 1 



N 



log(2) M 



AN, 



(34) 



where we used (128(1 . This shows how fragile the initial 
pure entangled states are when they contain a large num- 
ber of photons. Even when a tiny fraction e of photons is 
lost, the number of photons lost may still be appreciable 
if \a\ 2 is large, and this may lead to a large degradation 
of entanglement. 

When e increases the entanglement loss per photon 
decreases. Using (|25|) [since we assume e\a\ 2 is not small 
we cannot use (|28l) : we do assume í « 1 now] we plot 
úEn /dN as a function of log 10 e in Fig. 1. The 4 curves 
plotted in Fig. 1 are shifted versions of one another: in 
fact, they are really only functions of AiV = e| a| 2 . 

The approximation l|31|l . too, depends only on AN, 
not on S or e separately. In addition the change in entan- 
glement Ai?jv does not depend on M either since both 



Em{p) and Sjv(|$)($|) depend on M in the same way. 
To verify this behavior and to see whether this behavior 
extends to parameter regimes not covered by the con- 
ditions í < 1 and e <C 1 we plot in Fig. 2 the change 
in entanglement AEn as a function of AiV for various 
different parameters. For large M we use the approx- 
imation (|26|l . but for M = 20 we use the exact result 
l|22p. The plot confirms the predicted behavior: for 5 not 
too large the change in entanglement is only a function of 
AiV, but not of M, or e or 6 independently. For S = 0.2 
(and M = 20) there is only a small deviation from that 
universal behavior. 

Fig. 3 plots, in contrast, not the change in entangle- 
ment but the entanglement Ejq{p) left after decoherence. 
The parameter regime considered is very different as e is 
no longer assumed very small. We do assume { <C 1 
so that we can use the approximation 125|) . From the 
approximation l|31|l to that approximation we see that 
En{p) no longer depends on M, provided M is not too 
small. That behavior, too, is confirmed in Fig. 3, where 
curves for M = 20 (recali that whenever M < 50 the 
exact result (|22H is used) are seen to overlap with one for 
M = 20, 000. For smaller vàlues of M there is a clear 
dependence on M, as is illustrated for M = 2,3, 5. From 
Fig. 3 one notices that at reasonable (i.e. realistic) vàlues 
of e the amount of entanglement left is less than one unit, 
even if we started out with a state with log 2 M ps 14.3 
for M — 20, 000. However, that result is not optimized 
for the number of photons |a| 2 in the initial state. In 
particular, since we chose 5 to be small, \a\ 2 is roughly 
equal to M 2 . 

In order to see whether smaller numbers of photons 
may leave more entanglement after decoherence (since 
one expeets the decoherence to have a less strong in- 
fluence when there are fewer photons even if the initial 
pure state does not quite have the full log 2 M ebits of 



FIG. 3: Logarithmitic negativity En(p) as a function of 
e — 1 — r\. There are two sòlid curves almost completely 
overlapping that correspond to M = 20 and M — 20000, 
demonstrating that the amount of entangiement left after de- 
coherence is insensitive to M when M is not too small. The 
dashed curve corresponds to M = 5, the circles to M = 2, 
the diamonds to M = 3. 
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FIG. 4: Logarithmitic negativity En(p) as a function of \a\ 2 . 
There are three pairs of curves, each pair consisting of a sòlid 
curve, correpsonding to r\ = 0.7, and a dashed curve (for 
r\ = 0.49). The three pairs correspond to M = 20, 30, and 
40, with the màxima of the curves with larger M moving to 
larger |q| 2 . 



entangiement), we plot in Fig. 4 E^(p) as a function of 
\a\ 2 for various vàlues of M. Here we cannot use the 
approximation (|25[l and hence use the exact calculation, 
which, as mentioned before, does restrict the calculation 
to M < 50. The plot shows that the optimum amount of 
entangiement remaining after decoherence is almost inde- 
pendent of the value of M, thus confirming the findings 
from Fig. 3. For example, for 77 = 0.7 one is left with 
about 0.92 ebits of entangiement for the optimum vàl- 
ues of \a\ 2 almost irrespective of how much entangiement 
the initial pure state had. On the other hand, there is a 
strong dependence on \a\ 2 . Finally, the plot also shows 
that the decohered states become more robust against 



FIG. 5: Similar to Fig. 4: En(p) as a function of |a| 2 for 
r\ — 0.7 and M — 2,3, 4, 5, 6, 7, 8, . . . 10, where the màxima of 
the curves move to larger vàlues of \a\ 2 for larger M. 

photon absorption losses. For, even if one started out 
with a very large amount of (pure-state) entangiement, 
losing 30% of the photons (77 = 0.7) leaves one with at 
best around 0.9 ebits of entangiement. But a further loss 
of 30% (so that rj = 0.49) of the photons merely decreases 
the entangiement by a factor of 2. 

In Fig. 5 we plot the same function but for smallcr 
vàlues of M. As it turns out, if one fixes ï? = 0.7 the 
largest amount of entangiement one can be left with, 
starting out with multi-dimensional entangled coherent 
states, is in fact E max ps 0.95 achieved by taking M = 8 
and H 2 ps 2.86. 

IV. CONCLUSIONS 

We considered the decoherence of multi-dimensional 
entangled coherent states due to photon absorption 
losses. Those states are bipartite entangled states of two 
modes, and are described by two independent parame- 
ters. One is a, such that \a\ 2 gives the average num- 
ber of photons in each mode. The other is M, an inte- 
ger indicating the maximum possible amount of entan- 
giement between the two modes, log 2 M. The amount 
of photon absorption is governed by an additional pa- 
rameter, e, the fraction of photons absorbed from each 
mode. We calculated the amount of entangiement left af- 
ter decoherence as a function of a, e, and M. Onc useful 
quantity that can be constructed from the above param- 
eters is AiV = e\a\ 2 , the number of photons absorbed 
from each mode. Another useful quantity is S defined as 
S = exp(— | et | 2 1 exp(27rz/M) — 1| 2 ), which determines the 
overlap between two coherent states of amplitude a and 
phase difference 2tt/M. In terms of those quantities the 
conclusions are 

1. For large M the entangiement disappears at a rate 
given by 2/log(2) pa 2.89 ebits per photon as long 
as AN < 1. 
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2. For large M, small S and not too large vàlues of 
AiV the change in entanglement is a function of 
AiV alone. 

3. For large M and AN, and for small 6 it is the re- 
maining amount of entanglement that is a function 
of AN alone. 

4. In general there is an optimum number of photons 
|a| 2 for fixed e and M that maximizes the amount 
of entanglement left after decoherence. For that 
optimum, S is not very small (around 0.2). 



As an example of point 4, we found that almost one ebit 
of entanglement can be obtained with 30% photon ab- 
sorption by using a little less than 3 photons per mode 
for the initial pure state and M = 8. The corresponding 
value of (5 is about 0.19. That is the maximum amount 
possible, but almost the same amount of entanglement 
(at least 95% of the optimum) can be obtained by us- 
ing almost any value of M and the appropriate value for 
H 2 » M 2 /30. 
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